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A new type of soliton with controllable speed is constructed generalizing the theory of slow-light 
propagation to an integrable regime of nonlinear dynamics. The scheme would allow the quantum- 
information transfer between optical solitons and atomic media. 
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Solitons 0,0 SB are stable wave packets occurring in 
many areas of the physical world, from tsunamis, flocks of 
migrating arctic geese |5( to light pulses in optical fibres 
[Ella]. Their stability stems from the non-linear balance 
between dispersing and focusing processes. As a feature 
of non-linear dynamics, the speed of a soliton may depend 
on its amplitude. For example, the height of a solitary 
water wave determines its velocity However, once a 
soliton is launched there is usually no further control over 
its speed. Here we show how to generate optical solitons 
in atomic media that can be slowed down or accelerated 
at will. Ultimately, this method will allow the storage, 
retrieval and possibly the manipulation of the quantum 
information of solitons 0, Ell a, El in media. This idea 
extends the simplest scheme [lfj, |ll| for slow-light prop- 
agation 0, 0] to a genuinely non-linear regime. 

In the pioneering experimental demonstrations [lolllll 
ll^ l of slow light atomic media store the shape of light 
pulses in the spin states of atoms in what is known as 
dark states [l^J- Usually, such experiments operate in a 
regime where the atomic spins deviate only slightly from 
a default direction, a regime of linear spin waves 
Here we consider strong spin modulations where, as we 
show, the non-linear dynamics of light and atoms creates 
polarization solitons, see Fig.Q]for an illustration. Unlike 
other solitons known so far, their speed can be controlled 
after they have been launched, in precisely the same way 
as for slow light. In quantum-information applications 
0, such slow-light solitons could complement the use of 
quantum solitons in fibres 0, @ with the advantage of 
storing quantum information in media and complement 
methods for quantum memory with the advantages 
of non-linear dynamics, in particular the intrinsic sta- 
bility of solitons. Strong spin polarizations can be im- 
aged by illuminating the sample from the side with reso- 
nant light of uniform polarization. Moreover, they could 
be manipulated with light or magnetic fields, possibly 
a step from quantum-information storage to quantum- 
information processing in atomic vapours. 

Consider a cell filled with atomic vapour .11] or a Bose- 
Einstcin condensate illuminated with light propagat- 
ing in one direction. Only three atomic levels shall inter- 
act with the light, an excited state and two de gen erate 
ground states, for example two hyperfme levels |lfl 1 1 1| . 




FIG. 1: A slow-light soliton is a polarization twist. To visual- 
ize the polarization state of light we use the Stokes parameters 
si, S2 and S3 (with S3 in z direction) defining a vector at 
each point in space z and time t. For a given time t we plot 
a strip spanned by z and the polarization vector in arbitrary 
units. When time progresses, the polarization twist moves 
through the medium with velocity 0. The figure illustrates 
the spatial polarization profile at time t = in spatial units of 
the soliton length JHJ- The Stokes parameters [lTfl are calcu- 
lated from the solution Q for a slow-light soliton with linear 
initial polarization, Rabi frequency = 0.5MHz, spectral pa- 
rameter £ + it) = 10.0 x exp(i0.47r) MHz, Qo = $0 = and for 
a sharp line with f(A) = S(A). 



As a consequence of the conservation of angular momen- 
tum, the two circular polarizations of light couple to ei- 
ther one of the ground states. This type of interaction 
between three-level atoms and light is called A configu- 
ration 18]. In the presence of light of constant but arbi- 
trary polarization the atoms relax by spontaneous emis- 
sion to spin states pointing in the opposite direction of 
the optical polarization, so-called dark states 0] , a pro- 
cess known as optical pumping. The light spin-polarizes 
the atoms. Once the atoms are in dark states, they follow 
the polarization of light by themselves, without sponta- 
neous emission, as long as their excited states are sparsely 
populated 0] . Suppose that the polarization of the light 
has varied in time and has imprinted a spatial modula- 
tion in the atomic spin distribution. Any subsequent 
light incident with uniform polarization will attempt to 
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re-polarize the atoms, but, in turn, is re-polarized itself. 
The tail-end of the spin profile is reset to the default 
polarization and the front-end moves. The light thus 
carries the spin modulation along the medium, possibly 
modulating it further, with a speed that depends on the 
intensity. When the light is switched off, the spin profile 
stands still, as long as it is not eroded by diffusion, de- 
polarizing atomic collisions or the spin precession in stray 
magnetic fields. Illuminating the atoms will continue to 
move the spin wave until it reaches the boundary of the 
medium where the stored spin profile emerges as a po- 
larization profile of light. The pioneering demonstrations 
[lfl Hlj | of stopping light have used this method 0] , 
but with an important qualification [l4j : One of the po- 
larization modes of light dominates the other. In this 
regime, the spin-perturbation waves are small, with lin- 
ear dynamics, and they move without changing shape 
(li) . In the case of strong polarization modulations, the 
balance between the two competing processes involved, 
light polarizing the atoms and the atoms polarizing the 
light, may generate stable polarization structures, soli- 
tons USB- 

To develop a quantitative description of slow-light soli- 
tons, we assume that the light propagates in positive z 
direction with amplitudes that vary slowly in compari- 
son with the carrier frequency w. We describe the left 
and right circular polarization amplitudes in terms of 
the Rabi frequencies fl + (t,z) and z), respec- 

tively, i.e. by the transversal components of the electric 
field strength divided by H and multiplied by the dipole 
moment k of the atomic transitions the light interacts 
with. We distinguish atoms with different detunings A of 
the transition frequencies from resonance with the light, 
caused, for example, by the Doppler effect of their ther- 
mal motion. We describe the quantum state of each 
atom in terms of the density matrices pit, z, A) consider- 
ing only the ground states |±) and the excited state |e) 
that interact with the light. In the slowly varying en- 
velope approximation |2JJ] the light amplitudes obey the 
approximative Maxwell equations 



~dt 



dz 



+ OO 



{e\p\±)vdA, g 



(1) 

Here n(z) denotes the atom-number density profile of the 
medium and f(A) the relative distribution of the detuned 
atoms. £0 is the electric permeability of the vacuum. At 
this stage we neglect any atomic relaxation, but we deter- 
mine later the condition when this is justified. Without 
relaxation the atomic density matrices evolve according 
to the Liouville equation in the interaction picture 



dp 

at 

H = 



-i[H,p], 
-A|e)(e| 



-£(^l±)< e | + ^|e><± 
± 



(2) 



Despite the apparent complexity and non-linearity of the 
Maxwell-Liouville equations illll'll . they belong to the 
class of integrable systems [lg. The key element of an 
integrable system is a Lax pair 0>0,0] of matrices U and 
V that generates the equations of motion as the compati- 
bility condition dU/d(—dV/dT+[U, V] = for a variable 
complex spectral parameter .which allows the use of the 
Inverse Scattering Method [2J, IjJ Li • For the three- level 
Maxwell-Liouville system (11121) the detuning A serves as 
the spectral parameter, r means the retarded time t—z/c, 
and £ denotes the spatial coordinate z/c. One verifies 
that the Lax pair is [2l| 



U = -iH, 



p(A')is(A') dA' 



(3) 



Therefore, the theory of slow light is integrable 
which comes in useful for finding analytic solutions, and, 
more importantly, which allows the existence of stable 
non-linear waves, solitons. 

Various methods for finding soliton solutions start 
from a Lax pair of type We use a modification of 
the Zakharov-Shabat dressing method [22| for a variable 
background field f)(r) of constant polarization with the 
atoms being in the corresponding dark states. We obtain 
the single-soliton solution for the field f2± and the atoms 
p = ip a ipi as 



<p+ 
Q 

Ipe 



£ — ir\ tanh Q 



£ + if) 



= Qo 



= $n 
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(4) 



in the limit £ 2 + i] 2 \£l\ 2 where £, 77, Qo, <&o are real 
constants and f2(r) is a complex function. One verifies 
that the results (0} solve the reduced Maxwell equation 
JTJ and the approximate equation (A13) of Ref. [19] that 
describes the atomic dynamics. In the solution the in- 
cident field is left-circularly polarized. We can, however, 
apply the transformation 



= B 



f2_ 



(5) 



with constant unitary matrices B to describe an arbi- 
trary initial polarization state. The results |@J are gen- 
eralizations of some known solutions [2j| to slow-soliton 
propagation in the presence of inhomogeneous line broad- 
ening v(A), variable spatial density n(£) and, more im- 
portantly, time-dependent background O(r). 
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Slow-light solitons turn out to be pure polarization 
structures that, above a critical strength, perform com- 
plete polarization twists, see Fig. ^ Let us discuss their 
properties based on the analytic solution (J3J. A single 
soliton depends on four real parameters £, rj, Qq, $o and 
on one arbitrary complex function f2(r) for which we re- 
quire that \fl\ 2 < |£| 2 + 
® that 



We obtain from the solution 



\n+(t,z)\ 2 + \n-(t,z)\ 2 = \n(t- z /c)\< 



(6) 



which indicates that |f2| 2 describes the total intensity 
and that any incident intensity modulation propagates 
through the medium at the speed of light in vacuum. 
The parameter Qq characterizes the position and $o the 
phase of the polarization deformation, whereas £ and rj 
constitute the complex spectral parameter £ + irj whose 
argument 9 describes the maximal polarization devia- 
tion. Polarization twists occur when \9\ > n/2. Beyond 
the soliton the light returns to the incident polarization 
state, acquiring a geometric phase [24| of 29. The atoms 
are in pure quantum states with probability amplitudes 
ij) a (t,z,A). In regions of nearly uniform optical polar- 
ization the atoms are essentially in dark states, but their 
ability to follow a changing polarization of light is limited 
due to the detuning A. Outside the medium the solution 
(@J for ,g = describes the shape of the polarization pulse 
required to launch the soliton. In the medium, the soliton 
propagates with velocity 



mi 



C 2(£2+T/ 2 ) 



gu(A) dA 

(£_A)2 +??2 



2ff 
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for v/c <C 1 and in the second approximate expression 
for £ 2 + rj 2 much larger than a characteristic detuning 
A 2 . We obtain from the solution (J3J the length scale of 
the polarization profile, the soliton length 



r,\n\ 



V9 



(8) 



Typically for solitons 0j @, 0, the length is adjusted to 
the spectral parameter. Atypically U El 0; the speed 
hardly depends on any parameters, except for the exter- 
nal light intensity in exactly the same way as for tradi- 
tional slow light 0|, allowing control over the soliton 
after it has been launched. 

A critical issue for slow-light solitons is the potential 
threat of absorption, in contrast to fast self-induced- 
transparency solitons 0, E] that excite and de-excite 
the atoms on times much shorter than the relaxation- 
time scales of the atoms. Consider a slow-light soliton of 
length l s traversing the distance I in the medium. The 
principal relaxation mechanism is spontaneous emission 
from the excited state |e). The population of \e) is largest 
for zero detuning where \ip e \ 2 ~ 2(cy 'g)lj 2 v at the peak 
of the soliton. We use this case to estimate the absorp- 
tion. First, we express g defined in the Maxwell equations 



||TJ in terms of the Einstein A coefficient [23. We get 
g = 3/(167r) cAnX 2 , where A denotes the optical wave- 
length A = 2ttc/ui. In practice 0,0] the excited state 
decays into three ground-state levels, not only the |±) se- 
lected by the propagation direction of light. We estimate 
the fractional loss tjl as multiplied by the propagation 
time I jv times the maximal excited-state population, and 
obtain 



_ 32tt ;a 
VL ~^X 3 ~1 2 



(9) 



In the experiments 0, ^| nX 3 has been in the order of 
1 or larger. Losses are thus negligible for solitons much 
longer than the geometric mean of wavelength and dis- 
tance travelled, leaving enough room for soliton prop- 
agation. This requirement resembles the condition for 
traditional slow light 0] where the pulses are restricted 
in frequency space, which amounts to a minimal length 
scale in ordinary space. 

Having established that slow-light solitons are feasible, 
we identify their quantum properties. Fundamentally, 
the solitons are quantum excitations of the electromag- 
netic field interacting with the atoms of the medium, but 
they are very close to the classical limit (@J. Their quan- 
tum nature manifests itself in fluctuations around the 
classical amplitudes. A classical soliton is stable against 
fluctuations except when the fluctuations occur in its 
degrees of freedom, its parameters, in our case in any 
a G {£, 77, Qo, $0}, because such parameter fluctuations 
transform the soliton into another stable soliton. The 
Goldstone modes generated by infinitesimal parameter 
changes define the relevant quantum modes of the soli- 
ton Therefore, we represent the quantum fields Cl± 
and ip a of light and atoms as 



da 

4>a = + ^2 (a -a). 



(10) 



The quantum properties of the soliton modes depend on 
the parameterization of the classical soliton, chosen here 
(@J such that they are particularly simple. We deduce the 
commutation relations of the mode operators a from the 
quantum field theory of ligh t |25| in the slowly- varying 
envelope approximation [2fJ. 

Quantum field theories may start from a funda- 
mental quantum commutator and a classical mode 
decomposition with an appropriate scalar product. 
Here we describe the field by two components A+ 
of the vector potential in the Coulomb gauge 2a 
and SI units with the equal-time commutator 2£j 
[A±(t, zi), s$<7 dA±(t 1 z 2 )/dt] — ih8{z\ — z 2 ) where a de- 
notes the longitudinal cross section of the medium. In 
the slowly-varying envelope approximation 20] the Rabi 
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frequencies f2± are, apart from a prefactor n/h, the com- 
ponents of the electric field strength —dA±/dt that prop- 
agate with the carrier exp(— iwr). Therefore we define 
f2-t = k/(2H) exp(-iujT) (iujA± — dA±/dt) and get the 
commutation relation in retarded time 



^±(^0,^(72,0] 



9 

ft LO 



(11) 



Inspired by Lai's and Haus' quantum theory of self- 
induced-transparency solitons [9j we define the scalar 
product of the fluctuation modes in the style of Poisson 
brackets 



{a, /?} 



1 

87 



E 



da dp da 8(3 



dr. 



(12) 

One verifies by straightforward but lengthy calcula- 
tions that the solutions Q are designed such that all 
{a, 13} vanish, except {Q ,S,} = ~{£,,Qo} = {^o,??} = 
— {77, <&o} = 1- These orthonormality conditions serve 
to deduce the commutation relations of the fluctua- 
tion operators from the commutator We re-scale 
the quantum spectral parameter as £ + if] — (£0 + 
if) ) k 2 lu/ (16S()ti ac). Following the procedure by Lai and 
Haus [9| we obtain the Heisenberg commutation relations 



i = [Qo,£o] = [®o,fjo] 

= [ioM = [Qo,$o] 



[Qo,»3o] = [$o,6>]. (13) 



The components of the quantum spectral parameter, £0 
and fjo , thus represent the canonical momenta of two 
independent Heisenberg pairs with Qq and $0 as their 
partners. This implies that the position Qq of the soliton 
is complementary to £0 and the phase $0 of the polar- 
ization twist is complementary to its magnitude 7/0 • The 
quantum properties of these two modes encode the quan- 
tum information [t| of the soliton, regardless of whether 
it is moving or frozen. The overall amplitude Q plays the 
role of the control beam [ToL Hl| . When the soliton is 
stopped by switching off the incident light the quantum 
information is stored in the medium, since the atomic 
modes dip a /da do not vanish when f2(r) — * 0. Switching 
on the light will set the soliton in motion and will release 
the stored quantum information in the polarization of the 
light emerging at the boundary of the medium. 

Summary. — The theory of slow light 0, 0] is in- 
tegrable Q, y, m ■ A slow-light soliton is a polarization 
structure propagating with a speed that is proportional 
to the total intensity of the incident light. Entire solitons 
can be stopped and retrieved in atomic media with their 
quantum information stored in atomic coherences. 
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